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GIDROTEKHNICHESROYE  STROITEL'STVO  In  Russian  No  9.  1966  pp  15-20 


(Article  by  G.I.  Chilingarishvlli,  candidate  of  technical  sciences] 

[Text]  The  problen  of  the  temperature  conditions  of  blocks  of  dams  under 
construction,  being  concreted  consecutively,  allowing  for  the  effect  of  the 
foundation,  causes  considerable  mathematical  difficulties.  Several  solutions 
are  known  for  this  problem,  obtained  due  to  simplifying  assumptions  made. 

In  the  work  by  R.  Glover  [1],  a solution  is  given  to  the  one-dimensional 
problem,  using  the  analytic  method  for  blocks  of  equal  height  h,  placed  one 
after  the  other  at  regular  intervals.  To,  on  the  assumption  of  the  equality 
of  the  thermophysical  characteristics  of  the  concrete  and  the  rock  of  the 
foundation  and  the  assumption  of  other  simplifications.  The  heat  release 
function  is  taken  in  the  form  of  the  formula 


where  9o  and  m are  the  invariable  numbers  for  the  given  concrete. 


(1) 


R.  Carlson  [2]  solved  this  problem  with  the  same  prerequisites  by  the  method 
of  nets  according  to  a four-point  Implicit  system,  which  necessitates  the 
Joint  solution  of  s large  nu^er  of  equations  in  finite  differences.  The 
work  by  I.  Babushka  and  L.  Heyzdik  [3]  gives  the  results  of  solving  a two- 
dimensional  problem,  obtained  by  the  method  of  nets.  The  same  method  ms 
used  for  the  one-dimensional  problem  in  the  work  by  V.M.  Shteynberg,  I.  Ye. 
Prokopovich  and  I.V.  Gol'dfarb  [4],  for  a more  general  case,  when  the  thermo- 
physical characteristics  of  the  concrete  and  the  foundation  rock  are  not 
identical.  This  problem  is  most  widely  studied  in  the  recently  published 
works  of  Sh.  N.  Plyst  and  L.B.  Ssposhnikov  [5,  6].  The  use  of  a digital 
electronic  computer  made  it  possible  for  the  authors  to  examine  the  problem 
with  more  general  prerequisites.  Specifically,  they  used  the  finite  dif- 
ferences solution  of  one-dimensional  and  plane  problems  with  various  thermo- 
physical  characteristics  of  the  concrete  and  the  foundation  and  the  relation- 
ship of  the  heat  release  not  only  to  the  time,  but  also  to  the  temperature. 


Figure  1.  Graphs  of  Taaparatura  Distribution  in  tbs  Blocks  and  Foundation 
Dspsndlng  on  ths  Initial  Bast  Conditions  of  tha  Foundation 

l~initial  hasting  of  tha  foundation;  2—72  hours  after  placing  tha  block  of 
tha  second  tier  (144  hours  fron  tha  nonant  uhan  the  first  block  was  placed) 
3— sane  with  no  initial  heating  of  the  foundation. 
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In  all  the  above-mentioned  solutions,  obtained  by  the  method  of  nets,  equa- 
tions In  finite  differences  with  normal  precision  were  used.  Achieving 
the  assigned  precision  with  practical  calculations  often  requires  a sub- 
stantial refining  of  the  net.  Practical  considerations  often  require  Its 
considerable  enlargement,  and  therefore  a reduction  In  the  degree  of  error 
of  the  method  of  nets  becomes  quite  Important.  This  can  be  efficiently 
achieved  by  using  equations  In  finite  differences  with  Increased  accuracy. 

This  method  was  used  In  the  work  by  the  author  [7]  to  solve  the  unldlmenslonal 
problem  of  the  heating  of  blocks  placed  consecutively  outside  the  rock  area. 
Presented  below  Is  a solution  to  this  problem,  taking  into  consideration  the 
effect  of  the  foundation  with  different  thermophysical  characteristics  of  the 
concrete  and  the  foundation.  In  order  to  reduce  the  calculations,  we  shall 
restrict  outselves  to  a discussion  of  the  case  In  which  the  heat  release 
Is  an  arbitrary  function  of  time  6(t),  but  does  not  depend  on  the  temperature. 
We  shall  also  assutse  that  the  blocks  have  an  Identical  height,  h,  and  time  of 
curing  To.  We  will  begin  the  accounting  of  the  time  at  the  beginning,  from 
the  moment  of  pouring  each  new  tier  of  the  block,  and  the  upper  face  of  this 
tier  we  will  take  as  the  beginning  of  the  x axis  irith  Its  direction  downward 
along  the  vertical. 

The  thermal  process,  with  a source  of  heat  available.  Is  described,  as  we 
know,  by  a nonhomogeneous  differential  equation.  To  bring  It  to  a homo- 
geneous one,  we  Introduce  for  the  blocks  the  transformed  equation  function, 
related  to  the’ unknown  temperature  function  by  the  relationship 

(2) 

where  m Is  the  number  of  the  tier  of  the  block, n Is  the  number  of  the  time 
Interval  To  from  the  moment  of  pouring  the  given  block. 

Obviously,  we  have  the  following  equality: 

We  will  designate  by  A and  k the  coefficients  respectively  of  the  heat  and 
temperature  conductivity  of  the  concrete;  Ap  and  k o the  saM  coefficients 
of  the  rock  of  the  foundation.  Let  the  Initial  distributions  of  the 
temperatures  be  given:  In  the  blocks  — t^  OO,  In  the  foundation — t^^  ( }d 
and  the  temperature  of  the  upper  face  of  the  blocks  t (t) . Then  »e 
problem  for  the  transformed  temperature  function  of  the  block  of  the  1st 
tier  t',  (x.  t)  and  the  unknown  temperature  function  of  the  foundation  tQ(x>  '^) 
may  be  mathematically  represented  in  the  following  fora: 
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In  the  designation  of  the  transforned  tenperature  function  t^,  with  the 
transition  to  relationships  in  finite  differences,  for  convenience  of 
recording,  we  omit  the  upper  index  and  carry  the  lower  index  up,  placing  it 
in  brackets.  When  using  the  four-point  explicit  net  scheae,  the  equation 
in  finite  differences  for  the  block  of  the  n-tier  may  be  represented  in  the 
form  of  the  equality 


<!?  - «!a . /-^+  0 - *>«au  + (U) 

where  the  first  lower  index  with  t indicates  the  nusdier  of  the  pitch  along 
the  X axis,  and  the  second — the  number  of  the  pitch  along  the  t axis. 
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and  at  the  same  time  the  period  of  time  To  should  be  a multiple  of  the  pitch 
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Assuming  cr*l/6,  from  formula  (11)  we  obtain  the  equation  in  finite  differences 
of  increased  precision: 
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in  this  case,  from  formula  (12)  we  have: 


(16) 


Retaining  the  equation  (IS)  for  the  foundation  (with  the  substitution  of  the 
index  m for  0)  and  the  value  of  the  pitch  in  time  Ax,  the  pitch  along  the 
X axis  is  determined  by  the  following  formula: 
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Ue  record  the  conditions  (6)-(8)  in  the  finite  differences  In  the  following 
form: 


- . 'm.  . • 

(18) 

(19) 

(20) 

With  respect  to  the  conditions  at  the  border  of  the  two  media  (9)  and  (10), 
of  them,  (9)  may  be  used  directly.  With  respect  to  condition  (10),  we  note 
the  following:  with  the  transition  to  the  equations  in  finite  differences, 

the  even  curve  of  temperature  distribution  is  replaced  by  a broken  one,  with 
the  salient  points  at  the  nodes  of  the  net,  reflecting  the  changes  in  the 
heat  flow  with  certain  Juusps;  condition  (10)  expresses  the  equality  of  the 
heat  flows  entering  and  leaving  in  the  contact  node,  which  is  valid  when 
Ax^  0,  but  cannot  be  assumed  with  calculations  for  finite  differences. 


An  analysis  shows  that  the  values  of  the  temperatures  at  the  contact  of  the 
block  and  the  foundation  should  be  expressed  by  the  formulas: 
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where  C and  Y is  the  coefficient  of  heat  capacity  and  the  volume  weight  for 
the  concrete; 

0)  and  Yo are  the  same  values  for  the  rock  of  the  foundation. 

It  is  easy  to  see  that  formula  (22)  represents  the  condition  (9),  recorded 
in  finite  differences,  and  formula  (21)  in  the  particular  case  when  6*0  and 
0«1  are  turned  into  a formula  of  the  type  (15). 

From  the  moment  when  the  block  of  the  second  tier  is  poured,  the  transformed 
temperature  functions  for  the  blocks  will  be  connected  with  the  unknown 
functions  in  accordance  with  relationship  (2)  by  the  formulas 

<1  ( «.  ■«)  — tU (Of.  t)  + •(  (i);  ^24) 

(,(>•,  t)  + 4i(t).  ^25) 

The  transformed  temperatures  at  the  contact  between  the  blocks  in  the  finite 
differences  are  expressed  by  the  following  formulas: 
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Formula  (15)  la  used  for  the  inner  points  of  the  blocks,  and  for  the  contact 
between  the  block  of  the  first  tier  and  the  foundation,  formulas  (21)  and 
(22)  are  retained,  with  the  index  n replaced  by  the  index  2n. 


Sample  calculation.  A block  l.S  m high  with  a uniform  initial  temperature 
“ 10*C  is  placed  ont  he  foundation,  the  surface  of  which  is  heated  to 
20* C,  and  the  temnerature  at  the  depth  dies  down  to  12 *C  in  accordance  with 
the  given  law  t^0)(x).  The  block  is  cured  for  the  time,  To*72  hrs  at  a con- 
stant temperature  for  its  upper  surface,  equal  to  t4®^*10*C.  The  coeffi- 
cients are  assigned:  k*0.003S  m^/hr,  ko~2k,  Xo“2A  and  the  function  of  the 

adiabatic  heating  of  the  concrete  according  to  the  formula 
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where  0o-48.56*C,  a>2.293,  6-0.313. 


(28) 


From  formula  (23)  we  find  8-1.4142.  We. assume  for  the  concrete  block 
Ax-0.5  m and  from  formula  (12)  when  a-  we  find: 

i — % i;  *t: 


for  the  foundation  from  formula  (17)  we  obtain  Axq-0.7071  m. 

The  sequence  of  the  calculation  for  the  1st  and  2d  tiers  of  blocks  is  shown 
in  Table  1.  For  convenience  of  calculation,  a conditional  temperature  scale 
is  adopted  «rlth  zero  at  the  level  of  the  initial  temperature  of  thb  block, 
equal  in  our  example  to  10*C.  First  the  column  J-0  is  filled  in:  for  the 

points  of  the  block  of  the  1st  tier,  in  accordance  with  condition  (18), 
zero  is  recorded  (according  to  the  conditional  scale),  and  for  the  points  of 
the  foundation  the  initial  temperatures  are  recorded,  calculated  from  the 
given  function  t§  (x).  On  line  1 (i-0),  according  to  formula  (20),  the 
values  of  the  function  @i(t)  are  recorded,  with  the  reciprocal  sign,  pre- 
liminarily calculated  -in  a separate  table  for  each  pitch  of  At.  Then  the 
columns  are  filled  in  successively  from  top  to  bottom:  for  the  Intermediate 

points  (i-1,  2,  4,  5,  6,  ...) — from  formula  (15),  and  for  the  contact  points 
(1*3) — from  formulas  (21)  and  (22).  After  the  filling  in  of  the  column  for 
all  the  pitches  i,  the  values  of  A and  B of  this  column  are  calculated. 

These  values  are  used  to  calculate  the  temperatures  at  the  contact  point  of 
the  next  coloan.  As  the  result  of  all  these  calculations,  for  the  points  of 


'■  *1 


I 

the  foundation,  we  obtain  the  unknown  temperatures,  and  for  the  points  of 
the  block — the  values  of  the  transformed  temperature  function  t'^.  The 
transition  to  the  unknown  temperatures  for  the  points  of  the  block  Is  made 
In  the  next  column  of  Table  1,  according  to  the  formula  (2)  for  Tq~72  hours. 
Then  the  temperature  Is  calculated  In  the  concrete  masonry  and  the  foundation 
after  the  block  of  the  2d  tier  la  placed,  on  the  conditions  (2)-o 

(according  to  the  conditional  scale).  Line  16  (1-0)  and  column  J~0  for  the 
block  of  the  2d  tier  Is  filled  In  Just  as  for  the  block  of  the  1st  tier,  as 
the  result  of  the  Identity  of  the  Initial  and  limit  conditions  for  both 
blocks.  For  the  points  of  the  block  of  the  2d  tier  and  the  foundation  In 
the  column  J>0  (lines  25-28  and  31-39)  the  Initial  temperatures  for  the  points 
of  the  block  of  the  1st  tier  and  the  foundation  are  transferred,  respectively, 
from  the  column  of  the  latter  and  j>6.  Lines  20-22  and  24  are  also  filled 
In  preliminarily. 

Entered  as  the  latter  In  the  column  J-0  are  the  values  A",  A*  and  B',  which 
are  calculated  from  the  formulas  Indicated  respectively  In  lines  23,  29  and 
30.  The  next  columns  of  Table  1 are  filled  In  consecutively  from  top  to 
bottom,  by  calculating  the  temperatures  in  the  Intermediate  points  of  the 
blocks  and  the  foundation  from  formula  (15)  and  at  the  contact  points — from 
Che  formulas  given  In  lines  19,  25,  28  and  31.  In  the  last  column  the 
transition  Is  filled  In  to  the  unknown  temperatures  In  the  blocks  of  both 
tiers  with  the  aid  of  formulas  (24)  and  (25)  for  Che  last  time  Interval 
(J«6).  This  transition  Is  not  required  for  the  points  of  the  foundation, 
since  for  them  the  values  of  Che  unkno%m  temperatures  are  obtained  directly. 

Figure  1 shows  Che  distribution  curves  of  the  temperatures  In  the  two  tiers 
of  blocks  and  In  the  foundation  at  the  end  of  the  period  of  curing  the  block 
of  the  2d  tier  with  To**72  hours  with  the  pitches  of  the  net  Ax*0.25  m and 

hours.  Marked  by  circles  are  the  results  for  curve  2 with  the  enlarge- 
ment of  Che  Intervals  of  the  net  to  Ax*0.5  and  At*12  hours  (according  Co  the 
data  from  the  last  column) . 

The  curves  In  Figure  2 show  the  temperature  distribution  In  the  same  blocks 
with  various  Cheimophyslcal  characteristics  of  the  concrete  and  the  rock  of 
Che  foundation.  The  calculations  for  these  graphs  were  made  from  the  same 
net  with  pitches  Ax*0.2S  m and  At«3  hours.  Plotted  with  Che  same  pitches 
of  Che  net  In  Figure  3 were  the  temperature  dlstrlcuclon  curves  In  the  area 
near  Che  rock  of  a concrete  block  6 m high,  constructed  with  blocks  1,  1.5, 

2 and  3 m high.  These  curves  express  the  temperature  distribution  at  the 
end  of  the  curing  of  the  blocks  of  the  last  tier. 

From  the  results  of  the  calculations  made,  "envelope"  curves  may  be  plotted, 
for  each  variant,  of  the  greatest  temperatures  attained  at  each  pitch  for 
this  period  of  existence  of  the  block.  The  ordinates  of  these  "envelopes" 
for  the  areas  of  a block  1 m high,  adjacent  to  the  foundation,  for  the 
period  studied  by  us  (288  hours)  are  shown  In  Table  2. 
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Table  1 (continued) 
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Figure  2.  Graphs  of  Temperature  Distribution  Depending  on  the  Thermo- 
Physical  Characteristics  of  the  Concrete  and  the  Foundation 


For  values  k*0.0035  m^/hr;  1 — when  B“1.4142;  2 — when  B«l;  3 — when  B-0.7071; 

4 — for  k*0.0026  m2 /hr  and  B*0,707l 

The  accuracy  of  the  calculations  made  and  the  choice  of  the  optimal  pitches 
of  the  net  are  of  practical  interest.  The  evaluation  of  the  accuracy  of  the 
results  according  to  Runge's  principle  [8,  9],  obtained  for  the  above  graphs 
with  the  pitches  of  the  net  Ax*0.2S  m and  At-  3 hours  shows  that  the  degree 
of  error  everywhere  remains  less  than  0.1*C,  l.e.,  has  the  order  of  a round- 
off error,  usually  assumed  with  thermal  calculations  for  structures.  There- 
fore, the  results  with  these  pitches  of  the  net  can  practically  be  regarded 
as  fully  precise.  When  the  pitches  of  the  net  are  enlarged  to  A)^0.5  m and 
Atwi2  hours,  the  greatest  deviation  of  the  result  for  the  curves  given  In 


Figure  1 Is  0.2**C  (at  the  contact  between  the  blocks),  and  the  average 
deviation  remains  less  than  O.l'C.  Here  the  number  of  nodes  of  the  net 
resulting  from  Its  enlargement  is  reduced  eight  times  and  the  amount  of 
computation  work  decreases  accordingly. 
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Key: 

1.  Height  of  point  from  the  foundation,  In  m 

2.  Highest  temperatures  attained  ("C)  In  concreting  with  blocks 

of  a height  of 

In  practice  It  Is  often  sufficient  to  make  the  temperature  calculations  of 
the  .concrete  blocks  with  an  accuracy  up  to  0.5-l®C.  With  the  suggested 
method  this  accuracy  may  be  ensured  even  with  a more  substantial  enlarge- 
ment of  the  net  up  to  m an<i!lT*48  hours. 

We  note,  however,  that  to  avoid  greater  errors,  should  not  be  accepted 
over  1/3  the  height  of  the  block  in  the  Initial  period  of  heating  the 
masonry.  In  Figure  3 the  nodal  temperatures  for  the  fourth  variant  of 
constructing  the  block  (two  tiers  of  blocks  each  3 m high),  obtained  with 
a calculation  with  pitches A)fl  m and  At~48  hours,  are  marked  by  circles. 

The  intermediate  points  between  them  every  0.25  m,  connected  by  a dotted 
line,  are  obtained  by  Interpolation  at  Intermediate  Intervals  of  by 
Bessel's  formula,  which  Is  most  convenient  for  the  existing  Inflections  of 
the  curves,  and  In  the  uppermost  Interval — by  Newton's  formula:  It  Is 

Impossible  to  use  Bessel's  formula  here.  For  both  formulas  we  were  limited 
to  differences  of  the  second  order.  It  can  be  seen  from  Figure  3 that  the 
precise  curve  4,  In  practical  calculations  may  be  completely  replaced  by  a 
dotted  one,  and  Its  construction  by  the  above  combined  method  takes  only 
a few  hours  of  work  by  one  calculator.  With  this  substantial  reduction  In 
the  computation  work,  there  Is  no  longer  any  need  to  connect  them  with  the 
use  of  a digital  electronic  computer.  The  greatest  deviation  of  the  dotted 
curve  from  the  precise  one  here  Is  less  than  0.6*C  (at  a height  of  5 m from 
the  foundation  of  the  block),  and  the  average  deviation  Is  less  than  0.1*C. 
To  achieve  the  same  order  of  precision  on  the  basis  of  equations  in  finite 
differences  of  normal  precision  would  require  that  pitches  be  taken  at  least 
no  larger  than  Ax*0.25  m,  and  At*9  hours  and  the  number  of  nodes  of  the  net 
would  increase  over  20  times. 
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Figure  3.  Graphs  of  Taaparacura  Distribution  In  Concrete  Block  6 a High 

12  Days  After  the  Beginning  of  the  Concreting,  Depending  on  the 
Height  of  the  Blocks  h With  an  Identical  Rate  of  Building  the 
Structure  Ah"0.S  a per  Day 

For  values  k«k  "0.0035  a2/hr  and  B"l;  l--h"l  a;  2--h«1.5  a;  3 — h"2  a; 


The  facts  presented  above  Indicate  the  effectiveness  of  the  method  suggested 
for  determining  the  temperatures  in  dams  under  construction. 
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